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SECTION  I 


INTRODUCTION 


C -band  and  X-band  radars  are  commonly  used  to  track  satellite  and 
powe red -racket  vehicles.  The  resulting  measurements --range,  azimuth 
angle,  and  eievation  angle --are  processed  in  a  navigation  computer  to 
obtain  the  ephemeris  or  trajectory  of  the  target  vehicle.  In  general,  the 
radar  is  located  on  the  ground  and  the  target  vehicle  is  above  the  earth 
troposphere  («25  km).  Hence  the  radar  measurements  are  altered  from 
their  geometric  values  by  tropospheric  refraction.  The  actual  (geometric) 
value  for  the  range  may  be  7  to  600  ft  less  than  the  measured  range . 

Although  there  is  negligible  change  in  the  azimuth  angle,  the  actual 
(geometric)  value  of  the  elevation  angle  may  be  30  mrad  less  than  the  mea¬ 
sured  value.  Thus,  it  is  necessary  to  correct  the  radar  measurements  for 
tropospheric  refraction  before  they  are  processed  in  the  navigation  computer 

Simple  formulae  that  were  used  for  the  Atlas  and  the  Gemini  launch 
vehicle  radar  corrections  are  given  in  Refs.  1  and  2.  These  equations  are 
included  in  Appendix  B.  The  Atlas  real  time  correction^  is  in  error  by 
0.  01  deg  when  Eq  =  6  deg  and  by  05  deg  wnen  Eq  =  3  deg.  The  somewhat 
more  sophisticated  Gemii  i  correction^  reduces  the  error  to  0.  005  deg 
when  Eq  »  6  deg  and  0.  02  deg  when  Eq  »  3  deg.  In  both  cases,  the  error 
increases  rapidly  below  Eq  *  3  deg.  In  addition,  the  range  refraction  correc 
tion  error  is  often  r*100  ft  at  these  low  elevation  angles.  Clearly,  errors  of 
this  magnitude  can  lead  to  significant  navigation  and  ephemeris  errors.  One 
of  the  most  serious  errors  is  the  increasing  velocity  error  that  arises  from 
the  increasing  refraction  errors  that  result  when  simple  iormulae  are  used. 


The  ionospheric  refraction  correction  is  essentially  negligible  at  C -band 
and  X-band  radar  frequencies. 
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This  error  buildup  can  provide  velocity  magnitude  errors  of  several  fps  and 
pitch  velocity  errors  up  to  50  fps.  Hence,  it  is  desired  that  the  refraction 
correction  minimize  this  error  buildup. 

A  simple  refraction  correction  using  a  twelfth  order  polynomial  in  Eq 
is  derived  in  Ref.  3.  Even  in  this  case,  the  accuracy  below  Eq  =  5  deg 
deteriorates  rapidly  for  lack  of  higher  order  terms. 

More  accurate  refraction  correction  and  analysis  is  often  accomplished 
using  a  ray  tracing  program.  These  programs  divide  the  troposphere  bet¬ 
ween  the  radar  and  the  target  vehicle  into  many  (RslOOO)  layers,  each  having 
constant  refractivity  so  that  Snell's  law  of  refraction  may  be  used  with  each 
individual  layer.  They  require  many  mathematical  operations  for  each 
correction  and  therefore  may  not  be  usable  for1  real  time  applications.  Even 
with  ray  tracing  programs,  it  is  necessary  to  limit  the  bending  in  each  layer 
to  less  than  1  mrad  in  order  to  preserve  the  accuracy  of  the  technique. 

An  excellent  procedure  fox  correcting  radar  measurements  for  radar 
refraction  is  given  in  Ref.  4.  This  method  curve  fits  the  actual  index  of 
refraction  profile  using  locally  measured  data  (Rawinsonde)  and  integrates 
over  this  profile  to  obtain  AE  and  AR.  Unfortunately,  the  number  of 
computations  involved  limits  the  technique  to  use  in  postflight  reconstruction 
activities. 

Simple  refraction  correction  algorithms  suitable  for  use  in  Real  Time 
Applications  are  developed  in  this  report.  These  require  only  the  speclfica- 

V 

tion  of  the  surface  refractivity  Ng  at  the  radax  location.  Reference  5 

provides  a  detailed  discussion  of  the  rationale  for  using  N  in  this  manner. 

s 

This  simplification  leads  to  some  unavoidable  error  at  low  elevation  angles 
(below  2  deg)  because  of  the  variation  in  surface  moisture.  The  surface 

Nq  *  W  +  Dg  can  be  measured  using  local  atmospheric  pressure,  wet  bulb 
and  drylsulb  temperatures  as  shown  in  Appendix  A. 
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refractivity  Ng  is  the  sum  of  wet  and  dry  terms  (Dg  and  W  ),  and  these 
have  differing  profiles  with  altitude.  This  error  could  amount  to 
c  =  (^Eo)ft,0.05  deg,  which  is  considered  acceptable  since  most  radars  do 
not  function  well  at  these  low  elevation  angles  (**2  deg  or  less)  because  of 
ground  multipath  reflections.  A  correction  could  be  devised  that  used 
W  and  D  but  at  the  expense  of  model  complexity.  For  this  reason  this 
procedure  is  not  followed  herein. 

The  true  range  refraction  correction  is  dependent  upon  the  radar 
slant  range,  because  the  geometric  distance  traveled  for  the  actual  radar 
ray  exceeds  the  geometric  slant  range  to  the  target  vehicle.  This  excess 
distance  increases  somewhat  as  the  slant  range  increases.  However,  the 
variation  is  less  than  2  ft  whenever  the  elevation  angle  is  greater  than  2  de^. 
Hence  it  will  be  assumed  in  this  analysis  that  the  range  refraction  correc¬ 
tion  is  independent  of  the  radar  slant  range. 

The  improved  refraction  correction  developed  in  tills  report  assumes 

the  entire  troposphere  is  contained  in  a  single  layer  with  constant  N  .  The 

s 

height  of  the  layer  will  vary  to  agree  with  empirical  data  discussed  in 
Section  II.  The  model  development  is  discussed  in  the  following  sections: 

II.  Data  Discussion 

III.  Single  Layer  Bending  Model 

IV .  Slant  Range  Correction 

V.  Elevation  Angle  Correction 

VI.  Summary 

VII.  Validation 
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SECTION  II 


DATA  DISCUSSION 


The  data  for  developing  the  improved  refraction  correction  uses  the 

Centred  Radio  Propagation  Laboratory  (CRPL)  exponential  atmosphere  model 

given  in  Ref.  6.  Data  are  included  for  the  following  values  of  refractivity 

N  ,  measured  elevation  angle  E  in  mrad,  and  height  H  in  km  of  the  target 
s  o 

vehicle: 


N  =  200,  252.9,  289,  313,  344.5,  377.2,  404.8,  450 
s 

Eq  =  0,  1,  2,  4,  8,  15,  30,  65,  100,  200,  400,  900  mrad 

H  =  10,  20,  30,  50,  70,  90,  110,  225,  350,  475  km 

Note  from  Ref.  6  that  the  total  ray  bending  T  becomes  constant  and  the 
radar  slant  range  correction  AR  tends  to  become  constant  at  heights  above 
^25  km.  The  refraction  correction  developed  herein  will  only  apply  when 
the  target  vehicle  is  above  this  altitude. 

The  CRPL  data  is  included  in  Tables  B- 1  and  B-2,  giving  the  values 

of  t  and  AR  for  200  £  N  £450  and  0  £  E  £  900  mrad.  Table  B-2  includes  the 

s  o 

values  for  AR  when  the  target  vehicle  is  at  three  different  altitudes.  The  top, 
middle,  and  bottom  entries  are  for  target  altitudes  of  50,  475,  and  225  km, 
respectively.  Note  that  a  computation  of  AR  when  Eq  =  90  deg  is  also  inclu¬ 
ded  in  the  last  column  of  Table  B-2. 

Empirical  formulae  will  be  developed  in  this  report  to  essentially 
curve  fit  the  data  in  Tables  B-l  and  B-2.  Note  that  when  the  bending  t  is 
known,  a  simple  procedure  allows  the  computation  of  the  radar  elevation 
angle  refraction  correction  e.  The  physical  geometry  of  the  single-layer 
model  is  discussed  in  Section  III. 
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SECTION  III 


SINGLE  LAYER  BENDING  MODEL 


The  geometry  of  the  single  layer  bending  model  is  shown  in  Fig.  1. 

The  single  layer  of  refractivity  N  and  height  h  is  shown  as  a  shaded  area. 

s  , 

The  radar  at  Point  A  measures  an  elevation  angle  Eq.  The  ray  is  refracted 
through  the  total  bending  angle  t  on  departure  from  the  layer  at  Point  B. 

The  ray  continues  to  the  target  vehicle  at  Point  C. 

The  equation  for  Snell's  law  in  spherical  coordinates  gives  the  following 
equation; 


^  a  cos  E  =  (a  +  h)  cos  8. 

3  0  1 

where 

=  1  +  N  X  10’6 

8  S 

a  =  radius  from  earth  center  to  the  radar 
9j  =  elevation  angle  of  ray  after  departure  from  layer 
The  law  of  sines  for  triangle  OAB  is  as  follows: 

where  R  ^  is  the  length  of  the  line  AB. 

The  law  of  sines  for  triangle  OBC  is 


(1) 


a  +  h 


R  -  Rj 
sin  (5  -  $j) 


\  ~  # 
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where 

r  =  distance  from  center  of  earth  to  target  vehicle 

R  =  geometric  slant  range  from  radar  to  target  vehicle 

0  =  elevation  angle  of  ray  at  target  vehicle 

The  last  term  in  Eq.  (3)  is  not  an  equality  because  AB  and  AC  are  not 
colinear.  However,  the  result  is  an  excellent  approximation. 

The  law  of  sines  for  triangle  OAC  is 


(4) 


Equations  (1)  through  (4)  form  the  basis  for  the  single  layer  model. 

The  first  step  in  the  analysis  is  to  evaluate  the  thickness  of  the  single 
layer  such  that  the  proper  bending  T  results.  This  is  accomplished1  by 
solving  Eqs.  (1)  and  (2)  for  h  after  eliminating  0j.  The  result  is 

cos  E 

o 

cos  +  T 

where 


tan  0  j 


(5b) 


The  results  of  the  computation  for  h  are  included  in  Table  B-3.  Values 

for  h  in  nmi  are  included  for  each  case  (Eq,  Ng)  for  T  +  0.001  deg,  T,  and 

t  -  0.001  deg.  This  technique  permits  an  evaluation  of  the  sensitivity  of  r 

to  the  layer  thickness.  Note  that  when  E  <15  mrad  the  layer  thickness  must 

o 

be  accurate  to  about  20  to  50  ft.  As  E  increases,  the  tolerance  in  the 

o 
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layer  thickness  permits  errors  of  several  miles  without  causing  the  erroi  in 
r  to  exceed  0.001  deg. 

The  layer  thicknesses  along  with  the  tolerances  are  now  plotted  in 
Fig.  B-l.  Note  that  the  layer  thickness  is  a  linear  function  of  Ng  when  the 
elevation  angle  is  a  constant.  The  best  fit  for  h  is  given  in  nmi  by 

h  =  4-2  +  T555-  (62° -Ns)  <*> 

where  K  is  a  function  of  Eq.  The  key  variable  in  this  computation  is  Ng. 

That  is,  this  computation  does  not  depend  upon  the  height  above  sea  level  of 
the  radar.  Hence,  regardless  of  radar  height,  a  =  3444  nmi. 

Equation  (6)  is  then  solved  for  K  using  the  data  in  Table  B-3.  The 
results  are  given  in  Table  B-4.  Note  that,  as  before,  values  are  calculated 
for  t  +  0.001  deg,  t,  and  t  -  0.001  deg.  The  resulting  values  for  K  are 
plotted,  along  with  tolerances,  in  Fig.  B-2.  As  before,  when  Eq  £  15  mrad, 
the  tolerances  are  small;  as  Eq-*90  deg  the  tolerances  become  large.  The 
best  fit  for  K/1000  is  of  the  form 


K 

1000 


0.0083  -  (X.  0072  e"0,0221  E0 


where  E  is  in  mrad. 
o 

This  completes  the  empirical  curve -fit  process  for  the  evaluation  of 
t.  A  summary  of  the  equations  follows: 


K 


8.3  -  7.2  e 


-0.0221 


E 

o 


h  =  1.2  + 


T55o  <620  ‘  Ns* 


(cont. 
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(7) 


fj.  a  cos  E  =  (a  +  h)  cos  9. 
rs  o''  1 

(a  +  h)  cos  (0j  +  t)  =  a  cos  Eq 
a  =  3444 

These  equations  were  used  to  compute  t  xor  each  of  the  N  ,  Eq  data  cases  of 
Table  B-l.  The  results  are  given  in  Table  B-5.  The  maximum  error  is 
M).05  deg  when  Eq  =  0.  Other  errors  are  as  large  as  0.005  deg,  but  most 
errors  do  not  exceed  0.001  deg  for  Eq  - 1  deg  and  252.  9  ^  Ng  <377.  2. 
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SECTION  IV 


SLANT  RANGE  CORRECTION 


The  slant  range  correction  AR  is  based  on  an  observation  by  C.  Brown, 
of  the  General  Electric  Company,  Syracuse,  N.  Y.  Brown  noticed  that 
when  AR  was  large,  it  was  nearly  proportional  »•_  the  ray  bending  t;  also 
that  when  t  -  0,  AR-7  ft.  Hence,  the  expression 


was  computed  for  the  data  listed  in  Tables  B-l  and  B-2.  The  results  are 

given  in  Table  B-6  and  plotted  versus  N  in  Fig.  B-3,  It  was  noted  that  slopes 

s 

of  (AR  -  7)/t  tended  to  be  constants  when  Eq  was  constant.  Since  the  slope 
was  evaluated  as  1/25,  the  following  expression 


was  evaluated  (See  Table  B-7  for  the  results).  Note  that  the  values  of  b  are 
closely  grouped  for  constant  elevation  angles.  Central  values  for  the  groups 
are  listed  as  b.  Values  for  37.2  -  b  were  then  plotted  in  Fig.  B-4,  which 
shows  an  excellent  straightline  fit  on  log-log  graph  paper.  Hence,  b  was 
evaluated  as 


b 


37.  2  -  0.232  E 


0.604 

o 


Values  for  AR  in  Table  B-2,  applicable  to  a  target  altitude  of  475  km,  are 
used. 
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These  values  also  are  listed  as  b-calc  in  Table  B-7. 

Hence  the  resulting  expression  for  AR  in  feet  becomes 


AR 


7  +  T 


37.2 


N 

s 

25 


0.232 


e0.6°4) 


(8) 


where  t  and  E  are  in  mrad. 
o 

Values  for  AR  were  computed  using  the  Ref.  6  values  of  t  from  Table 
B-i  and  given  in  Table  B-8.  Results  are  also  given  in  Table  B-8  for  AR 
when  the  equations  of  Fig.  2  are  used.  Note  that  all  values  are  within  1  or 
2  ft  over  the  pertinent  values  of  Ng  and  Eq.  Also,  the  errors  in  the  com¬ 
puted  values  appear  to  be  random.  That  is,  no  significant  time  varying  bias 
exists  that  grows  as  the  radar  elevation  angle  decreases.  Thus,  velocity 
errors  resulting  from  this  type  of  error  buildup  should  be  nil. 


As  discussed  in  Section  I,  the  approximate  range  refraction 

correction  given  by  Eq.  (8)  is  independent  of  the  radar  slant  range.  As 

shown  in  Table  B-2,  the  approximation  is  accurate  to  ±2  ft  when  Eq^2  deg. 

The  error  can  be  as  large  as  10  to  20  ft  when  E  ^O.  Equation  (8)  can  be 

modified  to  include  the  geometric  range  correction  AR  by  modifying  the 

& 

curve  fit  procedure  given  in  this  section.  Then  the  term  (AR  -  7)/t  would 
be  replaced  by 


AR  -  7  -  kARg 

T 


where,  from  Fig.  1 
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AR  =  R.  +  BC  -  R 
g  1 

=  Rj  +  BC  -  jiRj  cos  e  +  BC  cos  (t  -  e)j 
=  2ri  sin2  |  +  2BC  sin2  ~~ - 
s  2Rj  sin2  |  +  2(R  -  Rj)  sin2 

and  k  is  derived  from  the  best  curve  fit.  This  compensation  for  AR^  is  not 
implemented  in  Eq.  (8)  because  the  expression  is  not  planned  for  use  below 
Eo  =  2  deg. 


SECTION  V 


ELEVATION  ANGLE  CORRECTION 


An  expression  for  the  radar  elevation  angle  refraction  correction  e 
will  now  be  derived.  It  is  assumed  that  the  measured  radar  slant  range  R* 
and  Eo  are  known.  The  procedure  for  this  calculation  using  Eqs.  (1)  through 
(4)  is  by  no  means  unique.  It  is  felt  that  the  method  included  herein  is  one  of 
the  simpler  methods. 

Equations  (7)  and  (8)  are  assumed  to  have  been  calculated  so  that  h,  9,  T 
and  AR  are  known.  The  procedure  will  use  Eqs.  (2),  (3),  and  (4)  to  compute 
Rj,  9,  and  finally  e.  The  quantity  Rj  can  be  computed  from  Eq.  (2)  as 

Ri =  ^rs-8in(T  +  ei  -E„>  »> 

o 

Equation  (3),  written  as 


R  =  R"  -  AR 


(ID 
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Equation  (4)  then  gives 


cos  (8  +  T  -  e)  =  ~  ain  (T  +  0  -  E  ) 

A,  O 


(12) 


which  can  be  solved  for  e . 

A  sensitivity  analysis  was  performed  to  evaluate  the  need  for  Eq.  (11). 
It  was  found  that  an  error  of  1  km  in  slant  range  R  results  in  an  error  in  *  of 
0.  002  mrad  or  less  when  the  vehicle  is  above  rn  altitude  of  20  km  or  greater. 
Therefore,  the  radar  measured  slant  range  R  can  be  used  for  R  in  Eqs.  (10) 
and  (12). 
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SECTION  VI 


SUMMARY 


A  block  diagram  of  the  single  layer  radar  refraction  correction  equations 

is  given  in  Fig.  2.  The  Eq  and  R  are  measured  by  the  radar  in  radians  and 

feet,  respectively,  while  N  is  measured  using  the  procedure  described  in 

s 

Appendix  A.  The  computational  procedure  shown  in  the  block  diagram  results 
in  the  elevation  angle  refraction  correction  in  radians  and  AR  in  feet. 

For  some  areas  it  may  be  desirable  to  "fine-tune"  the  constants  in  the 
block  diagram  from  their  nominal  values  of 


K1 

=  8.  3 

K6  ■ 

37.  2 

K2 

=  7.2 

K7  * 

-0.232 

K3 

=  22.  i 

It 

co 

0.  604 

K4 

=  1.2 

K9  * 

7 

K5 

=  620 

Kt<r 

25 

in  order  to  agree  with  ray-tracing  programs  using  climatology  data  for  a 
specific  area.  This  is  easily  accomplished  by  varying  the  ten  K  constants 
and  noting  the  effect  on  e  and  AR.  However,  no  such  fine  tuning  has  been 
found  necessary  for  the  VAFB  area. 
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fe 


(to  be  subtracted  ^  subtracted 

from  E  )  from  R  * 

o 


Fig.  2.  Block  Diagram  of  the  Single- Layer 
Refraction  Model  Equations 


SECTION  VII 

VALIDATION 


The  accuracy  of  the  single- layer  equations  to  compute  the  elevation 
angle  refraction  correction  c  has  been  tested  using  the  data  in  Ref.  6.  The 
results  are  given  in  Table  B-9.  Note  that  in  all  except  extreme  cases  the 
accuracy  is  good  to  **  0.  001  deg  (=S0.  02  mrad).  The  AR  computation  test 
results  for  the  same  cases  are  given  in  Table  B-8.  These  are  accurate 
within186'  2  ft  except  for  the  extreme  cases. 

The  equations  of  Fig,  2  with  the  nominal  constant  values  have  also  been 
tested  for  three  Atlas  F  flights  launched  from  VAFB,  and  an  Atlas  SLV  3  flight 
launched  from  ETR,  The  results  are  given  in  Table  B-10.  The  data  on  the 
correct  values  for  e  and  AR  were  supplied  by  C.  Brown  and  obtained  using  his 
ray  tracing  program.  Note  that  essentially  all  range  corrections  are  within 
±  5  ft  and  all  elevation  angle  corrections  are  within  ±  0.  001  deg.  These  four 
flights  span  an  N  from  293  to  360.  Hence,  it  is  felt  that  the  single  layer 
refraction  correction  model  with  the  equations  given  in  Fig,  2  is  sufficiently 
accurate  for  precision  guidance  and  tracking  purposes. 

This  testing  also  provides  a  validation  of  the  computational  procedure. 

All  computations  were  performed  using  a  Texas  Instruments  SR  52  which  has 
12  decimal  f*  40  bits  plus  sign)  accuracy.  Hence,  double  precision  arithmetic 
may  be  required  on  the  navigation  computer.  Floating  point  errors  were  only 
noted  in  computing  T  when  Ng  =  450  and  0  ^  Eq  <  15  mrad.  Special  logic  can  be 
included  to  prevent  floating  point  errors  in  this  case. 

The  computational  procedure  will  give  excellent  results  (Ac  2  0.  001  deg, 
AR  2  2  ft)  when  the  target  vehicle  is  above  an  altitude  of  25  km  (2  13.  5  nmi) 
for  Eq  <  100  mrad.  When  Eq  >  900  mrad,  this  altitude  may  be  reduced  by  50% 
to  12.  5  km  (6,  75  nmi). 
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APPENDIX  A 

COMPUTATION  OF  N 

8 

The  surface  refractivity  Ng  may  be  computed  using  T^,  T^.,  and  P, 

where 


Td  =  dry  bulb  temperature,  °F 


0 

Tw  =  wet  bulb  temperature,  F 

P  =  barometric  pressure,  mbar 

(1013.  25  mbar  =  760  mm  Hg  =  14.  7  Ib/in. 


Then  N  is  the  sum  of  dry  and  wet  terms: 
8 


where 


and 


with 


N  =  D  +  W 

8  8  S 


D 

8 


=  77.6 


P 

T 


e  =  the  actual  water  vapor  partial  pressure,  mbar 
T  =  the  dry  bulb  temperature,  °K 


(A-  1) 


(A-2) 


(A-3) 


T  =  |  (Td  -  32)  +  273.  15 


A-  1 


Then 


Eshbach  (Ref.  7)  shows  that  e  may  be  evaluated  as 


e  =  e  (x)  -  0.  000367  P  (T_  -  Tw) 
s  JJ  w 

where  e  (x)  is  the  saturated  water  vapor  pressure  at  the  wet  bulb  temperature 
s 

T.^  in  millibars. 

Reference  8  shows  that  an  excellent  approximation  to  the  saturated  water 
vapor  pressure  in  millibars  is  given  by 

e  (x)  =  6.  11  X  107*  5x/(237.  3  +  x) 
s 


Tw  +  1539\ 
157!  / 


where 


x 


32) 


In  some  applications,  T^y.  is  not  measured.  The  surface  refractivity 
can  be  computed  using  Eqs.  (A - 1 ),  (A-2),  and  (A-3)  when  relative  humidity 
RH  is  given  in  place  of  T^y.  In  this  case,  the  value  for  e  required  in  Eq.  (A-  3) 
is  computed  as 


e 


RH 

100 


e8(y) 


where 


y  = 


-  32) 


A-2 


Note  that  e  (y)  is  the  saturated  water  vapor  pressure  in  millibars  evaluated 
s 

at  the  dry  bulb  temperature. 

In  other  applications,  the  dew  point  temperature  in  degrees 

Fahrenheit  is  measured.  Then  the  partial  pressure  e  required  in  Eq.  (A-3) 
is  given  by 


• =  s(2) 


where 


z 


9  (TDEW  "  32) 


A-3 


APPENDIX  B 


SIMPLE  REFRACTION  CORRECTIONS 


Atlas  missile  real-time  correction  (Ref.  1): 


AE  =  N  X  10~^  cot  E  )l  -  7; — c' fr"'" )  *  in  rad 
s  o  )  C_o  am  E  1 

(  h  c 

N  X  10"6 

AR  *  -  cl'ain'E''  in  km 
E  o 


N 

CE  =  log  N  +SAN 
s 


0. 005577N 


AN  =  -7.  32  e 


where 

Eq  =  elevation  angle  measured  from  the  local  horizontal 

S  =  slant  range  from  radar  to  vehicle  measured  in 
kilometers 

Ng  =  surface  refractivity  in  n-units 
=  (u  -  1)  X  106 

s 

fig  =  surface  value  (at  radar)  of  the  index  of  refraction 
Gemini  real-time  correction  (Ref.  2): 


AE  =  N  X  10’6  cot  E  I  1  - 
s  o  | 

AR  ~  same  as  above 

a  =  earth  radius  =  6373  km 


C^S  sin  E 
E  o 


1  + 


2a  sin  E 


■  |  »  1 


in  radians 


B-l 


SURFACE  REFRACTIVITY,  N$ 


Fig,  B-3.  Plots  of  (AR  -  7  ft)/r  versus  E  and  N 

os 


B-4 


ELEVATION  ANGLE,  E  ,  in  mrads 


Table  B-i.  Total  Ray  Bending  Data,  t  in  mrad,  as  Presented  in  Ref.  6 
as  a  Function  of  Ns  and  EQ  when  the  Target  Vehicle  is  Above 
the  Troposphere  ^25  km). 
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When  only  one  E  value  is  given,  the  difference  between  these  is  negligible. 


Table  B-3.  Height,  h,  of  the  Single  Layer  of  Constant  Refractivity  Ns  in  nmi 
such  that  a  Total  Bending  of  t  +  0,  001  deg,  T,  and  T  -  0.  001  deg 
Results  where  T  is  Presented  in  Table  B-l. 
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The  three  values  in  each  box  correspond  to  T  +  0.001  deg,  T,  and  T  -  0.  001  deg  where  T  is  the  true 
value  of  total  bending  from  Table  B-l. 


200  28.  1796  27.948  27.  393  26.899  2S.  538  |  24. 652  (  21.821  16.  176  0.951 


252.  ''  I  27.099  26.827  26.18  25.671  24.668  23.  782  I  2 1 . 248  I  16.944  7.788 
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Results  are  included  for  Ng  =  200,  252.9,  289,  313,  344.5,  377.2,  404.8,  and 
50,  and  for  several  elevation  angles  in  rad  and  radar  ranges  in  km. 
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Table  B-9.  (Continued) 
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